Stationary axial plasma equilibrium^ 
magnetic polariza^io^l 



fiiwof the 




Robert W. .Johnson 

Alphawave : Research 
^AtiiS^k, GA, USA 

- ^~~~n@alphawaveresearch. com 



April 13, 2010 

rds^rMsma equilibrium; plasma magnetization; polarization force 
PACS: 28.52.Av; 52.25.Xz; 52.65.K 



Abstract 

The standard equilibrium equation for magnetized plasma is extended to ac- 
count for the magnetic polarization force. A factor of the pressure gradient arising 
from the magnetic decomposition of the Hall term survives the limit of vanishing 
magnetization and is canceled by a contribution from the magnetic polarization 
force, such that the proposed equation reduces to the standard form in the free 
current limit. Comparison of the solutions for an axially symmetric plasma col- 
umn indicates agreement with the free current limit for vanishing magnetization 
and variation in the field and current profiles as the degree of magnetization 
increases. 



1 Introduction 

The full e ffect of plasma magnetization remains poorly understood. From ear 



plorations (IGrad and Rubinl . 1 19581 ; iRose and Clarkl . Il96ll ; ISha franov . Il966t iR oth. 



Vlasovl. Il968t iHafll. 119721) t hroug h modern developments (IHazeltine and Meiss 



y ex- 



1967 



2003 



Hazeltine and Waelbroeckl . 120041 ) . we contend that an important effect has been ne- 



glected in the analysis of plasma eq uilibrium, the magn etic polarization force, despite its 



experimental applications in fusion (IHaves et ajJjT985b ma gnetic fluids 



2005 ) and materials science fjCantor and O'Reillvl. 



2003 
2007 



19991: iRinaldi and Zahnl. 120021: lEl-Dib and Moatimidl.l2002h. b iophysics (iQi et alllNovember 2001 



"Cozul 


:a et a 


Keil 


2007: 



Ata ka and Wakayamal . l2002l ; IWang and Waka vamal. 120 03: Tz irtzilakid.l2005l:lRamachandrana and Lesliel 



2002 



Zahn and Pioch 



Colli et all 120 03: Takari et al. 



2003 : lAsail. 12004 : Maki and Atakal . 120051 : Ma et alll2006l : IOno et al. 
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Jin et all 120081 ) . Its nonlinear nature is often cited w hen n e glect 



ing pla s ma magnetiz a tion i n the usual introductory equilibrium equation ( jChenl . 11984 ; 
Dendyl . Il993t IStaceyl . 120051 ). which treats the combined bound and free current densi- 
ties as a single entity; however, we find that the reciprocal relationship between M and 
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H leads to simplifications, rather than complications, in the treatment of stationary 
equilibrium in a strongly magnetized plasma. Note that it is a weak applied field H 
that produces the strongest magnetization M for a given pressure p. 

Following a brief review of conventional plasma equilibrium, we will examine the 
nonlinear model for plasma diamagnetism and its relation to the magnetic polariza- 
tion force. Next we consider the stationary equilibrium equation with the inclusion of 
magnetization to the net force balance. With restriction to an axially symmetric config- 
uration, we develop a model in the rest frame of the plasma for a magnetically confined 
column supporting an axial current . While there are several d ifferences with the usual 
model for solar coronal loops [see ( Tsypin and Galvaol . 2005 ) and references therein], 
most obviously the presence of gravity and the curvature of the column, we will neglect 
such effects in order to get at the root of the magne t ization problem. Application to 
axial experimental configurations (jCarter and Maggsl . 120091 ; iPalmer and Walkerl . 120091 ) 
is apparent. Performing a numerical evaluation, we compare the solutions of the mag- 
netically decomposed equation both with and without the magnetic polarization force 
with that obtained in the free current limit, where we find agreement only with its 
inclusion. 



2 Conventional plasma equilibrium 



There are many ways to view a collection of free electrons and ions immersed in a back- 
ground electromagnetic potential. The two most common approaches address either 
the kinetic equ ations or th e fluid equations, and collisional transport theory is used 
to relate them (jFitzpatricki . 120081 ) . The kinematically conserved quantities are mass, 
momentum, and energy, each related to an increasing velocity moment of the trans- 
port equation. The conventional plasma fluid equations rely on the statement that 
"all charges and currents are considered to be free" rather than utilize the formalism of 
macroscopic electrodynamics, and of primary importance is the influence of the Lorentz 
force on an electrically neutral medium J x B. In its simplest form neglecting gravity, 
the stationary equilibrium net force balance equation (also known as the momentum 
equation) for the single fluid model, 



Vp = J x B 



represents the balance of the Lorentz for ce due t o pla sma currents and the collision 
force identified as the pressure gradient (jWoodsl . 12004 ) in the absence of chang es to 
the momentum through Newton's second law. Much technology (ILao et all Il985l ) goes 
into the solution of that equation in various circumstances. The geometric form of 
Equation ([1]) indicates that the current and magnetic field vectors point along isobaric 
surfaces; that is, J • Vp = and B • Vp = for J x B ^ 0. The use of a scalar 
pressure rather th an a gyrotropic ten sor is often considered adequate for the analysis 



of plasma devices (jStacey et all 120061 ) . and the distinction will not be considered here. 



One commonly works in the rest frame of the plasma so that the fluid velocity Vy 
makes no appearance. 
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One's mathematical model should reflect the nature and scale of the physical system 
considered. The fluid model in this article is most appropriate for the description of 
macroscopic systems of sufficient density and extent that equilibrium thermodynamics 
is applicable. While cold plasmas are easily confined by material boundaries, the goal of 
hot plasma confinement is to keep the ions sitting still (on average) by forcing electrons 
through the gas while applying magnetic fields from external sources. A plasma is 
often considered a perfect conductor, in that the more it conducts, the more it wants 
to conduct, as the resistance decreases when its temperature increases. It should also 
be considered a perfect diamagnet, with a susceptibility on the order of unity at high 
pressures, as those free charges are given much latitude for gyration in response to an 
applied magnetic field. 

Our approach reflects that of classical field theory, combining macroscopic electro- 
dynamics with continuum fluid mechanics, and we summarize our notation and assump- 
tions: fully ionized hydrogenic plasma of species s £ {e, i} with total particle density 
n = n e + rii, total pressure p = nT = Y2 s n s^s, mass density p m = nm = 'Y^ is n s m s , 
free current density J j = ^2 s n s esV s , vanishing charge density p e = ^2 s n s e s = and 
mass velocity p m V^ = ^2 s n s m s V s = 0, and isotropic pressure = 2Wj = T s where 
T s <— k B T s . Our units are SI, but we express thermal energy T in eV. Note the total 
particle density is twice the neutral plasma density n = 2nn- Our treatment of magne - 
tiza tion most closely f ollows that by iGriffithsl (119891 ) : lHazeltine and Waelbroeck! (120041 ); 
and iKaufmann et al.l (120041 ). except that we will be fully decomposing the equilibrium 
equation in terms of H and M, a procedure not found in the existing plasma physics 
literature. 



3 Diamagnet ism, f3 limit, and the magnetic polar- 
ization force 

The fluid dipole moment per unit volume for each species is taken as M s = n s fi s for 
fx s = —(W^/Bg)H s , where the field B s felt by a single particle of species s within the 
unit of volume is the net field less the particle's own contribution B s / /xq = H+M— fj, s = 
H + Mjt + a s M s , where k ^ s and a s = (n s — l)/n s is a unitless factor. For a sufficiently 
dense plasma, n s ^> 1 such that a s — > 1 and B s — > B. The total magnetization of the 
neutral fluid is given by the net dipole density, which for p = p/fio and h = H/H may 
be written as the total number of particles times their average dipole moment, 



M = ^M s = -^(p s /|H + M| 2 )(H + M) , 

s s 

= -hM = -hp/ {H - M) = 2n (/x e + /LtJ/2 



(2) 
(3) 



and has (IJohnsonl . 120091 ) th e physical solution M/H = [1 - (1 - Ap/ H 2 ) l l 2 }/2 as 
the plasma is diamagnetic ( Marshall and Goldstein . 196lh . Ultimately, the proper 
treatment of mag netization requires t he us e of quantum theory, in particular as to 
account for spin (IHalzen and Martini . Il985l ). From the form of the solution for M 
one can immediately read a limit on the ratio of kinetic to free magnetic pressure, 
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Figure 1: Virtual trajectories (dashed) for the guiding center in a curved magnetic field 
(solid). The magnitude of /x should be considered fixed and its direction tangential 
to the guiding center trajectory. The path \h along the field line defines a contour of 
/x • H, hence the force //oV(/x • H) is everywhere perpendicular to d\n and the work 
done no J^ 1 d\n • V(/x • H) is zero. 



/3h = 2p/H 2 < 1/2 for a dense plasma. In terms of the net field B, we have 
P B = 2p/(H - Mf = f3 H /(l - M/H) 2 < 2, and the ratio M/H is limited to 1/2. 
In the dilute fluid limit n — > 1/rf such that a = (n — l)/n = 1 — l/2n — > 1/2, we 
find M s — > fi s so that the limits (3h M/H — > 1 when T; = T e . One must be careful 
to define the appropriate unit of volume for a magnetized plasma, which should be on 
an order no smaller than the cube of the ion gyroradius as the description of M in 
terms of a bound current density requires at least a single orbit within the volume unit; 
beyond that scale one must address by trajectory charges passing through the region. 
An infinitesimal dipole /j, s immersed in a field B s exper iences a force f q = V(/a c -B, s ) 



on the basis of the general theo r y of magnetized mater ial ( iGriffithsl . Il989t ICohen et al. 



2003t lHaus and Melcherl . 120091 ) . iFeynman et al.l ( 119631 ) gives a nice derivation in terms 



of virtual work. This force felt by the guiding center is here generalized to the dipole 
density by F# = V(M ■ B). One may argue that the fluid should not exert a force 
on itself, and so we will also consider F# = yU V(M • H), with either denoted by F M . 
We will not yet attempt a kinetic justification for the presence of this term among 
the moments of the Vlasov equation, suggesting that it should be instated at the level 
of the collisional Boltzmann equation, but rather will explore the consequences on 
the equilibrium equation of its inclusion, appealing to the mathematically well defined 
limiting process used in the reductions demonstrated below. The effect of the magnetic 
force f s on an individual dipole may be understood heuristically as the force of constraint 
keeping the guiding center on track with its magn etic field line in the thin flux tube 
picture (IFerriz-Mas and Nunezl . 120031 ; ISomovl . 120061 ) . and the corresponding torque t s = 
li s x B s keeps the guiding center aligned, as in Figure [TJ The established formalism does 
not care how fj, s is to appear, only that it has appeared, and the analogous effects for 
ferromagnetic materials should be familiar to anyone who has ever handled permanent 
magnets in a laboratory setting or otherwise. Thus, J x B + Fm is our generalization 
of the macroscopic force densities ( iMelcherl . 11981c iRosensweigl . 119821 ) given by Lorentz 
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and Kelvin, Flk = /^oJ x H + /xoM ■ VH, and by Korteweg and Helmholtz, Fkh = 
J x B — H • HV/i/2. One may argue with our form for the magnetic polarization force 
on two grounds, as for finite gyroradius the expression for the force per particle f s is 
approximate and as we have chosen to differentiate the dipole density. However, our 
choice leads to convenient simplifications and may be considered a first step towards a 
full accounting of magnetization effects on plasma equilibrium. 



4 Magnetized equilibrium equation 

The equilibrium net force balance equation we consider is given by 

Vp = 3 P x B + F M , (4) 

and the magnetic field is determined by V • B = and V x B x = fioJx, where the 
subscript X appearing in Ampere's law identifies the appropriate current source for 
each component of the magnetic field B//xo = H + M. (Equilibrium in the presence 
of gravity requires the additio n of the gravitational force F g = p m g to the RHS of t he 
force balance equation above ( Yoshikawal . 1981 ; Throumoulopoulos and Tasso . 200lh .) 



We require of our proposed equilibrium equation that it reduce to the standard form 
= /io(V x H) x H in the limit of vanishing magnetization, which we identify as the 
free current limit M <C H such that J& <C J/. The total current density is the sum of the 
external current density and the plasma current density J = J ex t + J p , where the plasma 
current is the sum of the free and bound currents J p = J/+J& = Vx (H— H eit )+VxM, 
giving us the magnetized equilibrium equation 

Vp = no [V x (H + M)] x (H + M) + F M , (5) 

and we will show that the magnetic polarization force cancels a contribution from 
the magnetic decomposition of the Hall term which survives the limit of vanishing 
magnetization p -C H 2 , noting that the curl of H ex t is zero within the plasma as well as 
that the bound current, defined as the curl of the magnetization, remains divergence- 
free regardless of the geometry and includes the effects of the pressure gradient driven 
diamagnetic current and the curvature and VB drift currents V x M = —V x (p/S 2 )B. 
Note that bound charges and bound currents are physically distinct entities; the former 
are governed by Gauss's law and associated with the presence of a binding energy 
between unlike charges while the latter are governed by the Maxwell-Ampere equation 
and associated with the presence of charges undergoing microscopic circulation. The 
gyromotion of plasma particles constitutes just such a microscopic circulation and is 
amenable to treatment as the magnetization of a fluid. 

Using the vector identity (V x B) x B = (B ■ V) B — V(B ■ B)/2 to decompose the 
Hall term into curvature and gradient contributions, we consider first the limit of the 
curvature term. Denoting lim^^ by =>- we find M/H =>• and 

(B/fi ■ V)B/// = [(1 - M/H)H ■ V] (1 - M/H)H (6) 
= (1- M/F) 2 (H- V)H- (1 - M/H)U[U-V(M/H)\ (7) 
=► (H-V)H. (8) 
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The gradient term may be written 

-\v(B/^) 2 = -\v{H - M) 2 = -\v{H 2 + M 2 - 2MH) (9) 

jL jL jL 

= — -VH 2 [l + (M/H) 2 ] + VMH , (10) 
2 

and the first term reduces to —HVH. The second term, however, survives the limit 

2VMH = [H-^-A^'^VH + HV^-^-Ap) 1 / 2 ] (11) 

= [l-{l-4p/H 2 fl 2 ]HVH + HVH -HV{H 2 -4p) 1 ' 2 (12) 

= [2 - (1 - 4p/H 2 ) 1/2 ] HVH - H(H 2 - 4p)~ 1/2 V{H 2 - 4p) /2 (13) 

= [2 - (1 - 4p/H 2 ) 1/2 ] HVH - (1 - 4p/H 2 y 1/2 (HVH - 2Vp) , (14) 

such that VMH =>- Vp. The magnetic polarization force Fm / Ho =>■ — Vp serves to 
cancel this term so that the decomposed equation reduces to the standard form Vp = 
(H • V)H — HVH in the free current limit. Using F B = — Vp, the curvature term is 
balanced by the gradient of the energy density 2p+ fj, (H -M) 2 /2 = p+fi (H 2 -M 2 )/2. 
Note that a magnetic polarization force density of n f = Yl s would only enforce the 
reduction in a region of constant density where Vp = nVT. 



5 Numerical evaluation 

Restricting consideration to an axially symmetric plasma column with d/d6 = d/dz = 
embedded in a constant external magnetic field H ex t — H® xt z, the free current within 
the plasma is here supposed to be purely axial J/ = Jf(r)z, giving rise to an polar 
magnetic field H/ = Hg(r)0. The free field H = H ex t + Hf satisfies V • B = 0, leaving 
us to evaluate the net force balance. Using the parameter f M e {0, 1} to indicate the 
absence or presence of the magnetic polarization force, some algebra yields a differential 
equation for Hg, where using F# 

|- [2(1 + 2f M )p + H 2 + iJ(tf 2 - 4p)V 2 ] = _I 0^ 2 [H + ( ^ _ 4^1/2] 2 > (15) 

and using F H we have 

|- [2p + (1 + 2f M )H 2 + (1 - 2/ M )iJ(iJ 2 - 4p) 1/2 ] = -i (§) 2 [if + (# 2 - 4p)V 2 ] 2 )(16) 

which are invariant under the transformation Hg — > —Hg corresponding to the two 
possible orientations for the axial current. The unmagnetized equilibrium equation is 
given by 

^l + He ^Jk = A H 2 , 

dr dr r 9 ' 

and our problem is defined as: given p(r), find Hg(r). 



(17) 
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Figure 2: Magnetized solutions with (□) and without (0) the magnetic polarization 
force Fb compared to the free current model (+) for p = 10 18 eV/m 3 . Curvature is 
neglected on the left and included on the right. The magnetized model without the 
magnetic polarization force does not agree with the free current limit. 



The density of a solar coronal loop is on the order of 10 14 ~ 10 18 /m 3 , and tem- 
peratures are found between leV and Ike V, with an polar fi e ld of a few mT and a 
column- aligned field up to several Tesla ( Tsypin and Galvaoi 2005 ). Parametrizing 
the pressure of a plasma column of meter radius by p(r) = p (l — r a ) with constant 
a > and using units of eV/m 3 , we consider the parabolic a = 2 profile with a central 
pressure po ranging from 10 18 to 10 22 in an external field of lOOmT, ensuring that the 
central (3° H < 1/2. We start by comparing the magnetized solutions both with (/m = 1) 
and without (J'm = 0) the magnetic polarization force to that obtained by the free 
current model Hg(r) = p^r 2 with po = 10 18 in Figure [2j Neglecting the curvature term 
(•V off) introduces a factor a/2 compared to its inclusion (-V on) on the magnitude of 
the solution profile. The solution for the magnetically decomposed equation without 
the magnetic polarization force in no way resembles the free current model's solution, 
which differs only slightly from the solution with the magnetic polarization force, and 
it diverges when the curvature term is included. 

Next we compare the solutions using F B and F H (including magnetic curvature) to 
the free current model for p between 10 21 and 10 22 in Figure [31 At moderate (3jj the 
solutions for Hg are virtually identical, and only at extremely high (3^ do the magnetic 
polarization force models become distinguishable, where F# is more similar to the free 
current model. The magnetic polarization force models differ in whether Jf should 
be more or less than the free current model's value. Finally, we compute the bound 
current Jj = Hx V(M/H) — (M/H)V x H using the F# model for the same range 
of central pressures, shown in Figure HI The polar bound current generally exceeds 
the axial bound current, except near the core of the plasma column. At the highest 
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r (m) r (m) r (m) r (m) 

Figure 3: Solutions for the free magnetic field (a)-(d) and the free current density (e)-(h) 
using either F B (□) or F H (0) compared to the free current model (+) for p between 
10 21 and 10 22 . The difference in the field between the models becomes apparent only 
at extremely high (3h- The magnetic polarization force models differ in whether the 
central current should increase or decrease from the free current model's value. 




r (m) r (m) r (m) r (m) 

Figure 4: Polar (□) and axial (0) bound current density solutions with F# for p 
between 10 21 and 10 22 . Near the (3h limit the bound axial current approaches ±20% 
the free axial current's value. 

f3° H considered, the axial bound current ranges over ±20% of the axial free current's 
value, indicating that magnetization effects become pronounced as one approaches the 
(3h limit. 
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6 Conclusions 



In conclusion, we have determined that our conjectured form for Fm provides just the 
right cancellation required to take the magnetically decomposed H + M equilibrium 
equation over to its usual form in the free current limit M <C H for either choice of 
magnetic polarization force model. The extra Vp found upon taking the limit of the 
Hall term J x B stymies attempts to proceed with the reduction of the force balance 
to the free current limit without its inclusion. By providing the necessary balance, the 
magnetic polarization force Fm allows for the extension of the standard equilibrium 
theory into the strongly magnetized regime as the Ph limit is approached. 

When compared to solutions of the free current model for an axially symmetric 
plasma column, we find considerable agreement between the magnetic polarization force 
models over a wide range of 0h- Only in the limit 2p/H 2 — > 1/2 do the models 
become distinguishable, suggesting that experimental discrimination will be exceedingly 
difficult. By isolating the free and bound currents, the magnetized equilibrium equation 
provides a more complete picture of what is happening inside the plasma. 

A model for plasma based on macroscopic electrodynamics is conceptually simple: 
Maxwell's theory of fields det ermines the potential from the source A(J) through the 
geometric relation d*dv4 = J (jRyderl . Il985l ; iDavisl . Il970[ ). and Newton's (or Einstein's) 
theory of kinematics determines the source from the potential J{A) through the Ohm's 
law equation, which depends as well on the momentum K(J, A). Ionized media is no 
less a material than any other mechanistic system, as its constituents must obey the 
physics dictated by their local environment and their interactions with other members. 
While the field theory is more or less straightforward, the influence of those particle 
interactions on the kinematics is daunting to compute, and therein lies the difficulty 
and beauty of plasma physics. 
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